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The question of determining the underlying Fermi surface (FS) that is gapped by superconductiv- 
ity (SC) is of central importance in strongly correlated systems, particularly in view of angle- resolved 
photoemission experiments. Here we explore various definitions of the FS in the superconducting 
state using the zero-energy Green's function, the excitation spectrum and the momentum distribu- 
tion. We examine (a) d-wave SC in high Tc cuprates, and (b) the s-wave superfiuid in the BCS-BEC 
crossover. In each case we show that the various definitions agree, to a large extent, but all of them 
violate the Luttinger count and do not enclose the total electron density. We discuss the important 
role of chemical potential renormalization and incoherent spectral weight in this violation. 



The Fermi surface (FS), the locus of gapless electronic 
excitations in k-space, is one of the central concepts in 
the theory of Fermi systems. In a Landau Fermi liq- 
uid at T=0, Luttinger Q defined the FS in terms of 
the single-particle Green's function G _1 (k, 0) = and 
showed that it encloses the same volume as in the non- 
interacting system, equal to the fermion density n. In 
many Fermi systems of interest the ground state has a 
broken symmetry. Here we study states with supercon- 
ducting (SC) long-range order, where there is no surface 
of gapless excitations, and ask the question: Is there any 
way to define at T=0 the "underlying Fermi surface" that 
got gapped out by superconductivity? 

^From a theoretical point of view, this question is of 
relevance to all superconductors irrespective of pairing 
symmetry or mechanism. The answers turn out to be of 
particular interest for strongly correlated superconduc- 
tors, where the surpring effects that we find are large 
enough to be measured experimentally. Angle-resolved 
photoemission spectroscopy (ARPES) Q has emerged 
as one of the most powerful probes of complex mate- 
rials and has been extensively used to determine the FS 
in strongly correlated systems, often from data in the 
SC state Hi. One of our goals is to understand ex- 
actly what a T=0 measurement can tell us about the 
FS. This is especially important in the cuprates where 
the the normal state must necessarily be studied at high 
temperatures and does not show sharp electronic excita- 
tions, expected in Fermi liquids, in contrast to the SC 
state which does show sharp Bogoliubov quasiparticles. 
Our results are also of interest for a completely different 
class of systems: strongly interacting Fermi atoms y| in 
the BCS-BEC crossover 0, Here too the question 
of an underlying Fermi surface is of direct experimental 
relevance Sj. 

In this paper, we first show that Luttinger's original 
argument |l| cannot be generalized to the SC state, and 
this violation is related to broken gauge invariance 0. 
We then explore various criteria for defining the "under- 
lying Fermi surface" in the T=0 SC state, using prop- 
erties of the single-particle Green's function G(k, tu) [lfj 



directly related to experimentally measurable quantities. 
We present results for the two systems described above: 
(a) the d-wave SC state of the high Tc cuprates which is 
dominated by strong Coulomb correlations, and (b) the 
s-wave superfiuid state in the BCS-BEC crossover regime 
of atomic Fermi gases with strong attractive interactions. 
We will show that the various definitions lead to FS con- 
tours which are not identical, but nevertheless agree with 
each other to a remarkable degree. All of them violate 
the Luttinger sum rule (area enclosed equal to fermion 
density) and we obtain a detailed understanding of this 
violation: its magnitude is related to the SC gap function 
and its sign to the topology of the FS. 

Fermi surface criteria: It is perhaps not appreci- 
ated that the question of the "underlying FS" in the SC 
state is non-trivial, because in BCS theory the answer 
appears to be simple. In the BCS state one can look 
at 10] G BCS (k,z) = (z + &)/{z 2 - El) and ask where 
G BCS (k, 0) = 0. The resulting surface coincides with 
£k = 0, the normal state FS on which the pairing insta- 
bility takes place. Thus it is tempting to use G(k, 0) = 
in a more general setting to define the SC state FS. 
This is analogous to Luttinger's definition except G(k, 0) 
changes sign through a zero in the SC, instead of apple in 
the normal case. However, it is important to note |9( that 
there is no analog of Luttinger's theorem for SCs. One 
can write the Luttinger- Ward functional in terms of the 
Nambu Green's function matrix G, and try to generalize 
Luttinger's proof pj. However TrG only constrains the 
difference (n^ — n±), which is trivially zero in our case 
|l2j. and not the sum 13]. This is related to the fact 
that spin Sl° tal is conserved in the SC state but number 
is not. Thus one cannot show in general that the surface 
flfT | G(k, 0) = in the SC state encloses n fermions. We 
will come back later to why the Luttinger count never- 
theless seems to work in BCS theory. 

We next turn to various alternative definitions of the 
FS. (i) ARPES measures 01 the one-particle spectral 
function A(k, us) — — ImG(k, ui+i0 + )/TT and thus one way 
to define the "underlying FS" is to look at A(k, u> = 0) 
to map out the locus of maximum ARPES intensity. We 
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also describe below a closely related minimum gap locus, 
also motivated by ARPES experiments 0, 13 • (ii) In the 
cold-atom experiments it is possible to measure the 
momentum distribution, and therefore we also discuss the 
(somewhat ad- hoc but well defined) criterion n(k) = 1/2 
to define a surface that separates states of high and low 
occupation probabilities, (iii) We show below that the 
quasiparticle excitation spectrum, even in the strongly 
correlated SC state, is given by E k = (££ + | Akl 2 ) 1 / 2 
where £k is the renormalized dispersion and Ak the gap 
function. We then look at the contour defined by £k = 
to define the "FS" . In addition to comparing the contours 
obtained by various definitions, we also discuss the extent 
to which these results differ from G(k, 0) = 0. 
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FIG. 1: (a): Momentum distribution n(k) along (0,0) — > 
— > (tt, 0) — > (0,0) for three different doping levels. 
Comparison of renormalized mean field theory with varia- 
tional Monte Carlo (MC) results of ref. 0: (b) n(k) and 
(c) nodal Z(k) as function of x. 

(a) High Tc Superconductors: We describe 
the strongly correlated d-wave superconducting ground 
state and low-lying excitations using a variational ap- 
proach 0, 0] to the large U Hubbard model H = 
— Ylrr' trr'ct^Cr'a + U J2 r 'M'M on a 2D square lattice. 
We choose [l6| the bare dispersion ek = — ^7k + ^'Ak with 
7k = 2 (cos k x + cos ky) and Ak = 4 cos k x cos k y . We work 
at an electron density n = 1 — x with hole doping x <C 1. 
Our variational ground state is |^o) = exp(— %S)V\BCS), 
where \BCS) is the BCS wavefunction with d x 2_ y 2 pair- 
ing, the projection operator V = Y\ T (1 — n r fn r jj elimi- 
nates all double-occupancy and finite t/U corrections are 
built in through exp(— iS) |l5|. Here we present the re- 
sults of a renormalized mean field theory (RMFT) using 
the Gutzwiller approximation [l?! ] which are in excellent 
agreement (see Fig.^b,c)) with those obtained using the 
variational Monte Carlo (MC) which treats projec- 
tion exactly. The RMFT approach has advantages over 
MC for our present investigation since we get much bet- 
ter k-resolution and we can study spectral functions. 

In the RMFT we minimize (H — fiN) to obtain 
self-consistency equations for the gap function Ak = 



A(cos k x — cos k y )/2 and Fock shift Xk 18]. These in turn 
determine the BCS factors i>£ = 1 - u 2 . = (1 - £k/Ek)/2. 
The renormalized dispersion £k = 5t£k — Xk — A* incor- 
porates bandwidth suppression by the Gutzwiller factor 
gt = 2x/(l + x), the Fock shift Xk, and the (Hartree 
shifted) chemical potential fi. E k = (££ + | Akl 2 ) 1 / 2 is 
the excitation energy °f the projected Bogoliubov 
quasiparticle (QP) state |k f) = exp(— iS)V^fk\ \BCS), 
where 7^ = (ukC kT - v k c_ k |). 

The spectral function is of the form A(k,u>) = 
A coh (k,w) + A inc (k,w). The coherent part A coh (k,w) = 
Z(k) [u^5(uj - Ek) + v^5(uj + Eu)] where the quasi- 
particle residue Z(k) = g t - 2g st J2k' e ^' v k'/ U ~ 
2 ^£kl]k' v k'/ u with 9t = 2x/(l + x) and g st = Ax/(1 + 
x) 2 . The coherent weight Z decreases monotonically with 
underdoping, vanishing as x — > 0, as seen in Fig. ^c), 
while the incoherent spectral weight j4 lnc increases with 
decreasing x as required by rigorous sum rules |l9j . 
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FIG. 2: (a,b,c): Various "Fermi surface Contours" as a func- 
tion of hole doping x. The contours plotted are (1) £k = 0, 
(2) n(k) = 1/2, and (3) the non-interacting Fermi surface 
ek = (J- ■ Note that the "minimum gap locus" (see text) is 
indistinguishable from £k = 0. (d): The fractional difference 
5n/n between the area enclosed by £ k = and n plotted as a 
function of hole doping x. 

Renormalized dispersion: The form of the excita- 
tion gap Ek = + | Akl 2 ) 1 / 2 suggests that we identify 
£k = as the "underlying FS" . We must emphasize, that 
despite the BCS-like form of Ek, the theory neither as- 
sumes nor implies the existence of sharp quasiparticles 
in the normal state with a dispersion £k- In Fig. |21we 
plot the £k = contours for various x along with others 
which will be discussed below. We see that the £k = 
contours are hole-like - closed around (ir, it) - for small 
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x, but electron-like - closed around (0, 0) - for x > 0.16. 
The precise x at which the "FS" topology changes is a 
sensitive function of the bare parameter 16] t' /t. 

Minimum gap locus: We plot in Fig.E|A(k, 0) which 
is the zero energy ARPES intensity. We can neglect 
the incoherent weight at ui = and write A(k, 0) ~ 
Z (k)r / (E£_ + r 2 ) , where the 5-function is broadened with 
a small T. We now follow a procedure developed in an- 
alyzing ARPES experiments Various cuts through 
k-space are taken perpendicular to £k = 0. On each of 
these cuts we determine the location of the maximum 
A(k, 0), which is also the same as minimum gap E^ (ig- 
noring the negligible k-dependence of Z(k) on this locus). 
The locus of min E^ defines the "minimum gap locus" . 
We see from Fig. that the curve £k = and the "mini- 
mum gap locus" , although not identical, are very similar 
at every doping. In fact this difference is not visible in 
Fig. |21 as it is less than the width of the lines used. 
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FIG. 3: (Color online) ARPES intensity maps at zero energy 
A(k,0) for x = 0.15 (left) and x = 0.35 (right). 

Momentum distribution: We calculate n(k) = 
(cl Ck CT ) and find the result n(k) = Z(k)u^. + n(k), where 
first term comes from coherent quasiparticles and the 
second term n(k) = (1 - x) 2 /2(l + x) + 0{t/U) is the 
incoherent contribution |20j. The evolution of rik with 
doping x is shown in Fig. ^a). We choose the contour 
fik = 1/2 |2JJ to (somewhat arbitrarily) separate states of 
high and low occupation and its variation with x is plot- 
ted in Fig.H We see that unlike £ k = 0, the n(k) = 1/2 
contour does not exhibit a change in topology and is 
'electron-like' down to very low doping. This qualita- 
tive difference arises because £k = and the "minimum 
gap locus" depend only on the coherent part of the spec- 
tral function, while n(k) is an energy-integrated quantity 
that includes incoherent spectral weight 23] . 

Luttinger Count: We see from Fig. that the vari- 
ous "FS" contours enclose areas different from n, which 
is the area enclosed by the non-interacting FS ek = fi ■ 
In Fig. |2{d) we plot the difference between the area en- 
closed by £k = and n as a function of hole doping x. We 
see that, in general, this difference is non-zero when the 
system exhibits SC long range order (0 < x < 0.4) [j^ . 
For x < 0.16, the £k = contour is hole-like and we find 
an area enclosed greater than n, while for 0.16 < x < 0.4, 
this contour is electron-like and the enclosed area is less 



than n. The x > 0.4 ground state is a normal Fermi 
liquid and the Luttinger sum rule is valid |24| . 

A simple way to understand the variations seen in 
Fig. El^d), which include both Gutzwiller and Hartree- 
Fock renormalizations, is not immediately obvious. How- 
ever, the BCS-BEC crossover analysis below will give us 
clear insight into both the (small) magnitude of the vi- 
olation observed here and the relation of its sign to the 
"FS" topology. 

Zeros of G: Next we compare the FS contours ob- 
tained above with the surface G(k, 0) = 0, though we 
note that the latter is not of direct experimental rele- 
vance. ^From the form of A(k, ui) obtained above, we 
get 

/+oo 
du,A inC (k,Lu)/LJ. (1) 
-oo 

It is clear that £k = corresponds only to the first term 
G coll (k, 0) = 0, and not Elto a zero of the full G. From 
the sum-rule constraints |l9j on A lnc one can show that 
the integral above is necessarily negative. It then follows 
that the zeros of G correspond to £k > 0. This implies 
that for small x, where £k = is hole-like, G = gives 
an even larger violation of Luttinger count. 

(b) BCS-BEC Crossover: The evolution of a Fermi 
gas from a BCS paired superfluid to a BEC of composite 
bosons has now been realized in the laboratory ja] us- 
ing a Feshbach resonance to tune the s-wave scattering 
length a s . The dimensionless coupling g — l/kfa s can 
be varied from large negative (BCS limit) to large posi- 
tive (BEC limit) values, with unitarity (g = 0) being the 
most strongly interacting point in the crossover. We use 
the T=0 crossover theory 0, 0] to gain further insight 
into the question of the "underlying Fermi surface" . 

The structure of the T=0 Green's function 0] m this 
case is exactly the same as eq. Q, with an excitation 
spectrum = (££ + A 2 ) 1 / 2 . In the well-controlled limit 
of large dimensionality, treated within dynamical mean 
field theory [25|. Z is close to unity and the incoherent 
spectral weight is small even at unitarity. Thus, to make 
our point in the simplest possible manner, we work with 
Leggett's variational ansatz 0,H^]. Even at this level, 
where incoherent weight vanishes, the implications of the 
various FS definitions are very interesting. 

It is then easy to show analytically |27| that all the 
definitions investigated above, yield the same surface in 
k-space for all values of the coupling g — l/kja s . This is 
given by the bare dispersion fi 2 P/2m = fi(g), where the 
chemical potential fi strongly renormalized from its non- 
interacting value. Even in the weak-coupling BCS limit, 
/i is less than the non- interacting et by an exponentially 
small amount of order A 2 /e/ and the Luttinger count is 
violated. This violation becomes increasingly severe with 
increasing g: as the gap increases, n(k) broadens and fi 
decreases; see ref. |26|. Eventually on the BEC side of 
unitarity (g ~ 1), fi goes negative and the surface £k = 
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shrinks to k = 0, beyond which one enters the Bose 
regime. To summarize: the "underlying FS" does not 
enclose the total number density n, its volume decreases 
monotonically with g and for g greater than a critical 
value it is zero! 

This analytical solution is modified quantitatively by 
correlation effects beyond the Leggett theory, but qual- 
itative effects like the decrease in fj, and broadening of 
n(k) with increasing g persist, as also seen in both nu- 
merical jl^l and experimental studies of n(k) ||. 

We now see how the renormalization of the chemical 
potential in the presence of a SC condensate directly 
leads to a violation of the Luttinger count. For not too 
large attraction, the violation has a relative size (A/ef) 2 , 
and a negative sign for a particle-like FS, i.e., the under- 
lying FS encloses a smaller area than the non-interacting 
FS. To see how the sign changes for a hole-like FS,we 
look at the BCS-BEC crossover in a lattice model such as 
the attractive Hubbard Hamiltonian 29] . It is straight- 
forward to show, using a particle-hole transformation, 
that sign of the /x-renormalization reverses going from 
a particle-like to a hole- like FS. Thus we find that for a 
hole-like FS, the underlying FS in the SC state encloses a 
larger area than the bare FS. These are exactly the effects 
seen in the strongly correlated d-wave SC in Fig.|2[d). 

Conclusions: We have analyzed various criteria for 
the "underlying FS" in the T=0 SC state. We have 
shown that a "FS" deduced from a SC state measurement 
necessarily violates the Luttinger sum rule and does not, 
in general, enclose n fermions. We have gained detailed 
insights into the magnitude and sign of the violation. 
Our results are of most interest for the high Tc cuprates, 
where they can be tested in ARPES experiments, pro- 
vided one can independently determine the electron den- 
sity. The existing ARPES data 0, on LaSrCuO 
show a small violation of the Luttinger count with a sign 
change, consistent with our results. We should empha- 
size that our theoretical considerations makes no state- 
ment about the finite temperature non-Fermi liquid nor- 
mal states. 
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